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20.  continued 

of  the  Polynomial  Terrain  Model's  characteristics  and  construction. 

The  next  section  contained  the  development  plan  identified  for  production 
implementation  of  the  polynomial  modeling  technique,  and  the  remaining 
sections  report  on  the  status  of  various  phases  of  this  development. 

The  results  showed  that  the  Polynomial  Matrix  method  is  the  most  promising 
of  the  various  digital  terrain  formats  (DPT) . 
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Background . 

The  feasibility  of  the  automation  of  the  stereocompilation  process 
was  first  demonstrated  by  the  Bausch  and  Lomb  Optical  Co.  In  1950  under 
a US  Army  Engineer  Research  and  Development  Contract.  The  first  practical 
equipment  was  produced  soon  after  by  Pickard  and  Bums,  Inc.,  also  under 
a Corps  of  Engineer  contract.  In  the  following  years.  Army  and  Air  Force 
supported  research  resulted  In  a succession  of  more  accurate  and  efficient 
equipment.  This  work  culminated  In  the  development  of  the  production  system 
entitled.  The  Universal  Automatic  ^p  £ompllatlon  Equipment,^  UNAMACE, 
under  a US  Army  Engineer  development  contract  In  1964.  The  Defense  Mapping 
Agency  Topographic  Center  (DMATC)  has  had  the  work-horse  UNAMACE  systems  In 
production  since  1965  and  currently  has  seven  systems  operating.  These 
automatic  stereocompilation  systems  produce  100  elevation  measurements 
per  second,  and  represent  the  most  prolific  production  source  of  digital 
terrain  elevation  data  from  stereo  photography  In  the  world.  The  UNAMACE 
systems  operate  In  a profiling  mode  In  order  to  facilitate  the  production 
of  orthophotographs,  another  of  Its  considerable  capabilities,  which  are 
required  Intermediate  products  In  the  defense  map  and  database  generation 
tasks. 


These  systems,  however,  are  not  without  limitations  based  upon  the 
current  state-of-the-technology  of  automatic  Image  correlation  In  stereo 
photography.  These  limitations  can  be  generally  characterized  as  "a  lack 
of  reliability;”^*^  and, more  specif ically,  result  in  the  following  three 
problems.  First,  the  automatic  stereocompilation  process  is  more  precisely 
an  automatic  image  correlation  of  two  small  areas  on  photographs  taken 
from  different  viewing  stations.  Camera  orientation  and  non-zero  ground 
slopes  combine  to  distort  the  geometry  In  the  two  pictures  considerably  - 
making  the  Image  correlation  much  more  difficult.  Correction  circuits  or 
algorithms  are  employed  to  minimize  these  effects  but  they  still  definitely  j 

contribute  to  the  total  system  errors.  Second,  these  Instruments  frequently  | 

get  unrecoverably  lost  due  to  poor  Image  quality,  cloud  cover,  water  j 

bodies,  and  other  reasons  resulting  in  a necessity  for  mantial  Intervention.  ! 

Third,  the  automatic  equipment  output  has  non-negllgable  high-frequency 

measurement  errors  that  must  be  "smoothed"  or  filtered  before  orthophotograph  i 

or  database  generation.  This  problem. is  due  not  only  to  system  measurement  ; 

errors  and  the  area  correlation  technique  employed,  but  also  to  the 
correlation  occurring  on  the  tree  tops  or  building  roofs.  Originally 

1 

these  problems  were  solved  for  contour  production  by  manual  extraction  and 

smoothing  from  the  UNAMACE  drop  line  output.  The  UNAMACE  drop  line  chart  | 

f 

was  an  orthographlcally  correct  plot  of  greyshades  keyed  to  change  | 

unambiguously  through  elevation  transitions  at  a specified  contour  Interval. 

While  some  research  has  been  directed  at  correcting  these  problems  at  I 

compilation  time,  the  majority  of  work  has  been  directed  toward  their  | 

i 

I 

1 

i 
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solution  In  a digital  computer,  post-processing  fashion.  This  choice 
of  emphasis  was  due  In  part  to  the  downward  spiraling  costs  of  time  on 
general  purpose  digital  computers  and  the  upward  spiraling  costs  for 
special  purpose  electronic  hardware  systems,  and  In  part  due  to  the 
already  recognized  need  for  digital  processing  of  the  elevation  data  for 
contours.  It  was  felt  that  data  verification,  editing,  and  smoothing 
could  be  most  effectively  performed  during,  or  as  a part  oi^  the  automated 
contouring  task. 

The  first  efforts  at  digital  computer  post-processing  of  the  digital 
terrain  data  were  performed  by  personnel  at  the  Army  Map  Service  (now  the 

! 

i 

Defense  Mapping  Agency  Topographic  Center)  which  demonstrated  the  feasibility  1 

I 

of  digital  contouring  and  contained  the  first  use  of  the  least-squares  j 

technique  for  smoothing  UHAMACE  data.^*^  | 

The  Electromagnetic  Compatibility  Analysis  Center  (ECAC)  was  the  first 
agency  to  augment  terrain  elevation  profile's  definition  with  stream  and 
ridge  Information,  and  the  first  to  use  digital  data  compression 
techniques  to  compact  elevation  data  grlds.^ 

In  1969,  the  US  Army  Engineer  Topographic  Laboratories  (USAETL) 

Initiated  research  and  development  efforts  to  Improve  the  efficiency  and 

quality  of  the  automatic  contours.^  This  effort  was  the  first  reported  I 

work  In  which  the  area  smoothing  technique  was  utilized  to  smooth 

UNAMACE  data,  and  the  first  to  consider  the  digital  mosaicking  of 

separately  compiled  elevation  data  sets.  During  the  period  of  1971  thru 

1974,  DMATC  and  USAETL  supported  contract  efforts  at  the  University  of 
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Virginia  directed  Coward  automatic  production  of  near-cartographic  quality 

contour  lines.  * This  research  contained  the  first  simultaneous 

application  of  the  least-squares  approximation  technique  with  functional 

terrain  modeling,  as  well  as  the  first  efforts  in  elevation  data  editing 

based  upon  separately  digitized  planimetrlc  and  hydrographic  data.  These 

development  contracts  resulted  in  terrain  modeling  and  automated  contouring 

software  entitled.  Contouring  via  Che  Surface  Averaging  £oncept,  CONSAC. 

The  basis  of  this  software's  terrain  modeling  algorithm  is  the  sequential 

determination  of  a matrix  of  numerous,  locally  valid  low  order  polynomials 

by  least-squares  approximation.  Low  order  polynomials  are  preferred  so 

that  solution  of  terrain  characteristics,  such  as,  contour  lines,  terrain 

slope,  or  new  grids  from  the  model  is  optimized;  and  because  empirical 

evidence  has  shown  that  higher  order  polynomials  do  not  improve  interpolation 
13 

accuracy.  The  necessity  for  a sequential  algorithm  is  due  to  Che 

requirement  for  efficient  modeling  over  arbitrarily  large  areas.  Reasonable 

core  storage  and  computer  processor  time  requirements  can  only  be  met  if 

the  modeling  technique  is  sequential  and  has  local  definition  based  upon 

processing  of  only  local  information.  The  requirement  for  a globally  valid 

and  conslstant  terrain  model  is  achieved  through  utilization  of  the 

12 

Weighting  Function  Interpolation  Technique  (WIT)  on  the  locally  valid 
polynomials  which  comprise  the  terrain  model.  This  technique  assures 
reasonable  interpolation  between  neighboring  local  polynomials  and 
assures  the  continuity  necessary  for  smooth  terrain  representations, 
either  contours,  profile  lines,  slopes,  or  reflectance  greyshades.  The 
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•Ignlflcance  of  Clioa*  arguMUta  waa  alao  given  great  emphasis  by  Leberl 
in  his  empirical  study  of  interpolation  schemes  for  digital  terrain 
data. 


Introduction. 


In  the  early  70' s the  automated  data  gathering  systems  employed 
in  related  areas  such  as  Landsat  and  projections  of  digital  database 
requirements  for  archival  storage,  weapon  systems  and  environmental 
studies  at  DMA  focussed  Increased  emphasis  on  the  need  for  digital  data 
compression  techniques.  Motivated  by  this  trend,  the  compaction  Inherent 
the  CONSAC  polynomial  terrain  model  was  first  Investigated  in  1973 
and  found  to  result  In  a compaction  ratio  of  30  to  1 over  a 75  square  mile 
test  area.^^  This  technique's  compression  of  digital  elevation  data 
grids  was  compared  to  linear  prediction  and  fourier  transform  methods  in 
seven,  thirty-five  square  kilometer  test  areas  by  Crombie,  et.  al.  in 
1974.^^  This  study  achieved  compaction  ratios  of  27  to  1 over  the  test 
areas  using  the  polynomial  modeling  technique;  and  for  comparable  modeling 
accuracy  this  ratio  was  over  four  times  better  than  the  best  compaction 
achieved  by  the  fourier  transform  method  and  was  over  seven  times  the 
best  ratio  achieved  by  the  linear  prediction  method.  These  were  the 
latest  research  results  when,  In  1975,  the  Defense  Mapping  Agency  shifted 
emphasis  and  funding  from  research  to  development  of  a near-term 
production  implementation  of  digital  data  compression  of  terrain  elevation 
Information.  This  report  details  the  status  of  an  ongoing  research  and 
development  effort  at  the  US  Army  Engineer  Topographic  Laboratories,,  which 
has  been  directed  toward  this  end  since  that  time.  The  first  section 
Includes  a discussion  of  the  Important  data  characteristics,  its  major 
applications,  and  the  needs  for  compression,  as  well  as  those  areas  where 
compression  Is  not  warranted.  The  second  section  discusses  the  various 
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published  terrain  representations,  their  capabilities  and  limitations 
with  specific  emphasis  in  the  areas  of  application  compatibility  and 
data  compression.  The  third  section  presents  an  overview  of  the  Polynomial 
Terrain  Hodel's  characteristics  and  construction.  The  next  section 
contains  the  development  plan  identified  for  production  Implementation 
of  the  polynomial  modeling  technique.  The  remaining  sections  report  on 
the  status  of  various  phases  of  this  development. 


I 


I 
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Automatic  Stereocompilation  Data  Characteristics  and  Applications. 


As  mentioned  above,  the  most  significant  limitation  of  this  type  of 
data  Is  a lack  of  reliability.  If  taken  at  face  value,  this  statement 
can  be  very  misleading.  Of  the  three  factors  mentioned  in  the  background 
section,  namely  photograph  geometry,  system  losses,  and  system  measurement 
errors,  only  the  latter  affects  the  majority  of  the  data,  and  these  effects 
can  be  minimized  by  rigorous  mathematical  treatment.  Errors  due  to 
photographic  geometry  differences  and  system  losses  are  more  difficult 
problems,  but  occur  much  less  frequently,  usually  only  3 to  5 percent  of 
the  time.  Currently,  these  problems  are  handled  In  two  ways;  by  UNAMACE 
hardware  and  by  post  processing  on  a general  purpose  digital  computer. 

The  UNAMACE  allows  for  both  the  specification  of  adverse  areas  by  the 
operator  (in  which  the  machine  just  holds  elevation),  and  for  intervention 
and  retrace  of  operator  specified  profiles.  The  post-processing  is 
currently  a simple  surface  slope  filtering  algorithm  which  requires  user 
input  of  a maximum  allowable  slope  over  the  area  considered. Combination 

c 

of  these  two  capabilities  reduces  the  number  of  gross  errors  significantly 
but  does  not  completely  eliminate  them.  Small  Interactive  systems  are 
being  considered  which  will  permit  cost  effective  manual  editing  of  the 
remaining  problem  areas. 

The  characteristics  of  the  automatic  stereocompllatlon  elevation  data, 
which  are  Important  to  the  modeling  techniques,  can  be  summarized  as: 

(1)  is  obtained  on  a very  dense  grid,  typically  over  5,000  points  per 
square  inch,  and 
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(2)  has  a non^na^llgibla  randoaly  distributed  "noise"  level  due  to 


syataa  errors. 

It  Is  laiwrtant  to  note  that,  currently,  the  major  requirement 

dependent  upon  this  data  la  the  production  of  orthophotographs.  The  ^ 

i 

density  of  the  elevation  data  Is  directly  a result  of  this  application;  j 

J 

the  production  of  high  accuracy  and  high  resolution  orthophotographs  j 

depends  on  a dense  grid  of  continuously  varying  (e.g.  smooth)  elevation  ! 

values.  Another  important  application  dependent  upon  this  data  is  the 

production  of  contour  lines.  At  first  glance  this  application  would  seem 

to  require  a considerably  less  dense  grid  than  is  currently  employed. 

13 

Leberl*  s empirical  ptudy  seems  to  support  this  view  since  he  showed  that 
for  most  cases  elevation  accuracy  is  linearly  related  to  grid  spacing.  For 
example,  most  1:50,000  map  steets  are  contoured  at  a ten  meter  contour 
interval  and  Leberl 's  results  indicate  that  a grid  spacing  of  50  meters 
would  be  adequate  to  insure  the  one-half  contour  interval  accuracy  required 
for  class  A maps  in  the  most  rugged  terrain  studied.  Also,  numerous 
other  study  results  quoted  by  Leberl  indicate  that  the  currently  employed 
12.5  meter  and  25  Mter  spacings  are  far  too  dense  for  the  accuracy 
required  in  the  contours.  However,  as  was  mentioned  earlier,  the  automatic 
stereocompilation  data  has  randomly  distributed  non-negligable  "noise" 
even  after  the  blunders  are  removed.  In  order  for  the  least-squares 
algorithm  used  to  accurately  determine  the  elevations  for  contouring  in 
the  presence  of  this  noise,  as  many  points  as  reasonably  possible  must 
be  used.  This  is  true  because  the  expected  error  at  any  given  grid  point 
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is  decreased  as  the  number  of  neighboring  points  used  in  its  determination 
is  increased.  For  these  reasons,  the  optimal  contour  determination 
occurs  when  all  of  the  dense  grid  of  elevation  values  are  utilized.  This 
would  seem  to  run  contrary  to  Leberl’s  empirical  evidence  only  until  this 
and  the  fact  that  his  study  assumes  zero  error  in  the  grid  values  are 
considered. 

The  other  applications  Impacted  by  the  automatic  stereocompilation 
data  are  standard  formats  for  data  exchange,  weapon's  systems  and  archival 
storage.  The  currently  established  DMA  standard  format^^  for  exchange 
of  digital  terrain  elevation  data  Includes  the  uncompacted  grid  written 
in  a profile  format.  This  mode  was  determined  optimum  because  it 
minimizes  the  impact  of  partial  data  loss  in  transit,  it  minimizes  the 
complexity  of  software  required  by  different  user  communities,  and  It 
maximizes  the  data's  compatibility  with  existing  user  hardware  and  software 
application  systems.  Most  current  weapon's  support  systems  depend  upon 
elevation  data  grids  and  the  costs  associated  with  modification  of 
existing  fielded  systems  would  be  prohibitive.  Formats  for  archival 
storage  are  still  being  formulated  and,  in  fact,  this  study  is  directed 
toward  their  optimal  definition. 

From  initial  analyses  it  appears  that  existing  production  systems  for 
orhtophotographs  and  contour  lines,  data  exchange  formats^  and  current 
weapons  systems  are  operating  optimally  using  elevation  data  grids.  The 

t 

hardware  and  software  developments  necessary,  as  well  as  the  additional 
time  required  to  pack  and  then  unpack  the  data^  more  than  outweigh  any 
possible  benefits  of  having  fewer  bits  to  carry  through  these  systems. 
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These  arguments  will  continue  to  be  true  for  the  high  throughput 

♦ 

I production  system  of  the  future.  However,  the  future  weapons  and  tactical 

■ support  systems  could  substanclally  benefit  from  digital  data  compression 

j'  ' If  It  were  compatible  with  the  archival  storage  formats  of  the  day  and 

I 

"designed  In"  from  the  beginning.  Possible  benefits  from  compression 
Include  smaller  unit  size  and  Increased  ground  coverage  and  accuracies. 

One  often  quoted  and  Incorrect  "fact  of  life"  Is  that  use  of  compression 
technique  must  slow  down  the  response  of  any  given  system.  This  Is 
! generally  true  In  systems  that  have  negllgable  data  access  times  or  do 

i;  not  allow  for  some  "look  ahead"  or  prediction  with  pre-access  and/or 

preprocessing  of  the  data.  However,  data  access  times  can  be  Improved 

by  compaction  for  applications  with  relatively  long  access  times  by  I 

! 

i 

allowing  for  the  use  of  smaller  but  faster  storage  systems  If  the 
compaction  regeneration  Is  sufficiently  efficient.  This  point  Is  made 
; clearer  when  one  considers  that  on  the  CDC-6400,  for  example,  over  j 

150  different  16-term  bicubic  polynomials  can  be  accessed  from  a set  of  i 

over  6,000  core  resident  polynomials,  and  each  evaluated  at  different 
and  arbitrary  X,  Y values  during  only  the  average  disk  head  movement  time 
required  of  a single  disk  access.  Further,  (as  will  be  discussed  later)  { 

' 6,000  polynomials  represents  a very  large  area.  For  systems  that  allow  | 

prediction  of  the  next  set  of  data  need,  the  use  of  advance  access  j 

techniques  along  with  parallel  or  pipeline  microprocessors  provides 
for  the  additional  time  and  processing  power  necessary  to  regenerate 
compact  data.  The  utilization  of  compression  techniques  would  most 

I 

I 

a 

? 
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beneficially  Impact  the  archival  storage  of  digital  data.  This 
application  Is  not  as  sensitive  tn  compression/regeneration  time  and 
Is  quite  favorably  impacted  by  reduced  storage  size  requirements  and 
costs. 

The  environmental  applications  of  this  data  have  not  been  covered 
because  they  are  much  less  well  defined  at  this  point  In  time. 


Digital  Terrain  Fonaats  (DTF) . 


For  the  purposes  of  this  study,  a Digital  Terrain  Format  will  be 
defined  as  any  consistent  digital  representation  of  terrain  elevation 
values  over  a given  area.  A digital  representation  Includes  not  only 
the  actual  binary  Information  stored  In  the  computer,  but  also  the 
algorithms  necessary  to  reference,  unpack,  and  continuously  define  the 
data  for  use.  Note  that  using  this  definition,  the  digital  terrain 
model  (DTH)  or  equally-spaced  elevation  data  grid  Is  not  a DTF.  This 
digital  storage  strategy  has  an  Implied  access  algorithm  but  only 
becomes  a DTF  when  a method  for  consistent,  continuous  definition  Is 
also  specified,  such  as  linear  or  spline  interpolation. 

This  section  will  present  and  discuss  the  currently  proposed  digital 
terrain  formats  with  special  emphasis  on  their  application  compatibility 
and  compression.  The  digital  terrain  format's  characteristics  summarized 
In  table  1 will  be  utilized  to  motivate  this  discussion. 

The  elevation  matrix  Is  Included  In  the  table  for-  comparison  and,  as 
discussed  previously.  Is  the  currently  utilized  format.  The  binary 
storage  format  In  this  case  is  Just  the  elevation  values.  The  storage 
overhead  Is  minimal  because  the  horizontal  location  of  each  elevation  Is 
Implied  by  Its  position  within  the  grid;  only  the  position  and  orientation 

of  one  point  In  the  grid  and  the  spacing  are  needed  to  use  the  Information. 

13 

The  Interpolation  algorithm  used  Is  arbitrary,  as  proven  by  Leberl,  since 
the  Interpolation  accuracies  of  the  various  schemes  are  essentially  equal. 
The  application  compatibility  has  been  discussed  previously  for  this  DTF, 


13 


and  the  low  relative  implementation' cost  and  zero  compaction  ratio 
are  obvious. 

The  Coded  Elevation  Matrix  uses  a Huffman  coding  of  the  Differential 
Pulse  Code  Modulation  (DPCM)  technique  applied  directly  to  the  elevation 
matrix  profiles. This  approach  has  been  applied  to  digital  image  data^^ 
and  Is  being  Investigated  at  USAETL  for  application  to  elevation  data. 

This  Investigation  represents  an  application  of  existing  well-knom 
errorless  digital  data  compression  techniques  to  elevation  data  grids. 

The  other  characteristics  of  this  format  are  essentially  the  aame  as 
for  the  elevation  matrix  except  for  the  eight  to  ten  compaction  achieved 
(based  upon  preliminary  results  of  reference  18) . The  application 
compatibility  is  only  given  a good  rating  because  of  the  unpacking 
required  to  utilize  the  data. 

The  next  three  formats  considered  are 'commonly  based  upon  the  use 

of  Irregularly  spaced  data,  a quite  old  and  logical  approach  for  data 

gathering  In  a manual  mode,  especially  If  the  data  extrema  are  the  points 

of  greatest  Interest,  and  will  be  discussed  together.  The  first  of  these, 

termed  Irregularly  spaced  data-A,  Is  currently  being  Investigated  by 
19 

Feucker.  This  approach  involves  the  greatest  binary  storage  overhead 
because  It  not  only  requires  three  values  per  elevation  point  (x,  y,  z), 
but  also  requires  up  to  six  pointers  or  links  per  elevation  point  In  a 
secondary  data  structure  needed  for  efficient  access  and  use  of  the 

I 

elevation  Information.  In  this  format  the  overhead  storage  required  to 
use  the  date  exceeds  the  storage  needed  for  the  data  by  a factor  of  two. 
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Little  work  has  been  done  todate  on  the  effect  of  different  Interpolation 
algorithms  for  this  type  of  data,  but  It  is  suspected  that  results  similar 
to  those  found  for  grids  will  result.  The  Irregular  spaced  data-B  refers 
to  work  by  Whltten^^  which  Involves  the  use  of  orthogonal  polynomials.  T'.ie 
irregularly  spaced  data-C  format  refers  to  Hardy 's^^  multiquadric 
equations  of  topography,  which  defines  the  elevation  surface  by  determining 
a function  composed  of  a sum  of  quadric  surface  equations  (the  polynomial 
equations  for  paraboloids,  hyperboloids,  etc).  In  these  formats  (B  and  C) 
only  functional  coefficients  are  stored  and  the  overhead  Is  minimal 
since  the  structure  of  the  modeling  function  is  Implied. 

Of  significance  are  the  facts  chat  the  orthogonal  polynomial  coefficients 
can  be  determined  analytically  (closed  form  solution  equations)  while  Che 
mulClquadric  equations  of  topography  requires  numerical  solution  (inversion 
of  a matrix  the  size  of  the  number  of  elevation  points) . Considerable 
computational  savings  could  be  realized  by  the  multlquadrlc  approach  if 
the  elevations  were  constrained  to  lie  on  an  orthonormal  grid,  then  the 
bilinear  form 

H - XAY^ 

and  Its  efficient  solution  algorithm 

-1  “1  T 

A - X ^H(Y  y 

21 

as  presented  by  Schut  could  be  extended  and  utilized.  The  interpolation 
algorithms  for  the  B andc  formats  are  defined  by  their  respective  functional 


! 
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approach.  The  application  compatibility  of  these  three  formats  were 


I 

all  judged  poor  because  of  the  poor  efficiency  of  grid  generation,  poor 
area  to  area  continuity,  and  the  difficulty  encountered  in  producing 
this  data  from  the  current  elevation  data  sources.  Hardy  does  address 
the  continuity  problem,  but  his  solution  of  redoing  adjoining  areas  as 
a new  larger  model  is  not  feasible  for  a global  representation.  The  poor 
efficiency  of  grid  generation  for  formats  B and  C is  due  to  the  evaluation 
time  for  the  extremely  long  polynomials  that  result  when  areas  of  reasonable 
(mapsheet)  size  are  considered.  The  compaction  ratios  achievable  by 
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Irregularly  spaced  data  are  based  upon  Mark's  and  Peucker's  works,  which 
represent  the  only  comparison  of  Irregularly  spaced  points  and  grids 
available.  Mark's  work  falls  short  of  an  actual  empirical  comparison. 


and  Instead  makes  use  of  the  assumption  that  grid  errors  are  linearly 
related  to  grid  spacing  to  get  his  results.  Further,  Mark's  crltera  for 
comparison  were  not  elevation  values  but  were  the  mean  slope,  terrain 
relief,  and  hypsometric  integral  over  the  test  areas.  The  relationship 


between  these  measures  and  the  one  most  desired  here,  elevation  accuracy, 
is  not  clear,  however,  they  represent  the  only  data  available  at  this  time. 
Mark  found  ratios  of  Irregular  to  grldded  storage  space  of  47.2  to  1,  2.3 
to  1,  and  1.4  to  1 using  Peucker's  approach.  Peucker  notes  quite  correctly 
that  the  highest  ratio  (due  to  terrain  relelf)  is  the  most  unrealistic 
since  measurement  of  extreme  value  differences  is  not  an  appropriate 
measure  of  grid  precision.  The  other  two  ratios  can  be  translated  to 
irregular  to  grldded  points  ratios  by  removing  the  overhead  of  Peucker's 
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and  this  results  In  ratios  of  14  to  1 and  8.4  to  1.  These  ratios  Bust 


not  be  weighted  too  heavily  since  they  do  represent  projections  based  on 
comparative  crltera  not  suitable  for  this  study.  Peucker's  one  test 
comparison  resulted  in  a 1.3  to  1 ratio  when  the  overhead  Is  Included.  The 
Implementation  costs  are  all  high  due  to  the  complexity  of  the  basic 
database  generation  required  by  these  formats,  and  their  poor  source  and 
application  compatibilities.  Peucker's  technique  received  the  highest 
relative  Implementation  cost  rating  due  to  the  additional  complexity 
required  in  the  binary  storage  overhead  definition,  e.g.  finding  all  the 
links. 

The  last  entry  In  table  1 Is  the  previously  discussed  Polynomial  Matrix 
format  which  defines  a sparce  grid  of  low  order,  locally  valid  polynomials 
as  the  basic  digital  storage  format.  As  In  all  matrix  or  grid-  approaches, 
the  binary  storage  overhead  Is  negllgable  because  the  horizontal  position 
of  each  data  element  Is  Implied  by  Its  position  In  the  matrix.  The 
Interpolation  algorithm  used  Is  the  Weighting  Function- interpolation 
Technique  (WIT) , which  utilizes  weighting  functions  to  produce  a globally 
valid,  smooth,  and  continuous  terrain  model.  At  any  given  point  the 
global  model  Is  defined  by  a locally  valid,  low  order  polynomial  which  can 
be  efficiently  solved  for  a smooth  data  grid,  contours,  slope,  or  other 
desired  terrain  characteristic.  The  accuracy  of  this  approach  has  been 
demonstrated  on  an  extensive  test  area  (over  600  square  kilometers)  In 
Cache,  Oklahoma.  The  source  compatibility  Is  good  because  the  polynomial 
matrix  Is  sequentially  computed  directly  from  profile  subsets  of  the 
stereocoBpllatlon  output  grid  using  the  efficient  least-squares  algorithm. 


I 
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The  compaction  ratio  listed  is  based. upon  new  results  that  are  reported 
in  later  sections  of  this  paper,  but  which  basically  involve  the  use  of 
standard  errorless  digital  data  compression  techniques  on  the  coefficients 
of  the  local  polynomials.  The  30  to  1 compaction  ratio  reported  for  this 
technique  earlier  in  this  paper  can  be  explained  as  follows;  a grid  of 
polynomials  is  computed,  which  is  roughly  90  times  as  sparce  as  the 
original  data  grid.  Since  each  polynomial  has  an  average  of  three  coefficients, 
this  results  in  the  30  to  1 ratio.  Further,  digital  compression  of  the 
coefficients  results  in  at  least  an  80  to  1 ratio.  The  moderate  implementation 
cost  is  due  to  the  relatively  high  source  and  application  compatibility. 

Grids  can  be  very  economically  generated  because  the  DTF  is  a simple  low 
order  polynomial  locally. 


Polynomial  Matrix  Format. 

The  polynomial  matrix  format  consists  of  a sparce  grid  of  locally 
determined  and  locally  valid  node  polynomials.  These  polynomials  are 
called  node  because  there  Is  one  polynomial  centered  over  each  node 
In  the  sparce  grid.  The  processes  Involved  In  their  generation  and  use 
(as  well  as  this  format's  Important  characteristics)  will  be  discussed 
In  this  section  In  chronological  order. 

The  output  of  the  automatic  stereocompllatlon  equipment  Is  a computer 
magnetic  tape  containing  numerous  contiguous  DTM's,  each  In  sequential 
profile  format.  Batch  processing  software  Is  utilized  on  a large  scale 
digital  computer  to  mosalck,  transform,  and  regrld  these  DTM's  Into  a 
single  consistent  DTM  which  covers  the  area  of  Interest  In  the  desired  map 

O'! 

projection  or  cartesian  co-ordinate  system.  For  a 1:50,000  mapsheet  this  i 

results  In  a DTM  composed  of  approximately  1,000  profiles  containing  2,200  ! 

points  each,  with  a ground  spacing  of  25  meters  between  the  profiles  and 

12.5  meters  between  points  In  the  profiles.  Further  processing  may  be 

utilized  to  edit  the  DTM  based  upon  separately  digitized  planlmetrlc  detail, 

12 

such  as  lake  boundaries,  rivers,  drainage,  and  rldgellnes.  It  Is  at  this 
point  that  the  polynomial  modeling  process  begins,  currently  based  upon  the 
batch  processing  software,  COMSAC  II,  on  a large  scale  digital  computer,  e.g. 
a CDC-6400  or  UNIVAC- 1 108 . Profile  subsets  with  approximately  fifty 
percent  overlap  of  the  DTM  are  sequentially  read  Into  core  memory  and  then  , 
square  ground  areas  are  sequentially  defined  along  each  profile  subset,  also 
with  approximately  fifty  percent  overlap.  (Currently,  these  overlapping 
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square  ground  areas  contain  rectangular  data  sets  of  17  by  33  of  the  DIM 

points  representing  0.16  square  kilometers  on  the  ground  for  1:50,000 

mapsheet  areas).  Each  square  ground  area  then  receives  the  following 

processing.  (Note  that  since  each  subarea  receives  Identical,  Independent 

processing  this  algorithm  contains  a high  degree  of  parallelism  that 

could  be  taken  advantage  of  on  one  of  the  parallel  processor  digital 

26 

computers,  such  as  the  Goodyear  STARAN) . First,  the  subarea  Is  approximated 

with  a "node"  polynomial  using  the  least-squares  crltera.  The  weighting 

employed  In  this  approximation  Is  as  specified  by  the  Maximum-Error 

Local  Approximation  Theorem. This  theorem  specifies  the 

most  liberal  (maximum)  error  distribution  allowed  In  local  or  "node" 

approximation  so  that  a specified  approximation  tolerance  can  be  achieved 

when  the  WIT  algorithm  is  utilized  to  form  the  final,  smooth,  and  continuous 

approximation.  At  this  point,  surface  slope  filtering  and  simple  statistical 

measures  are  utilized  to  locate  extreme  errors  In  the  data  and  to  measure 

the  goodness  of  fit.  The  effects  of  the  bad  data  on  the  approximation  are 

«. 

efficiently  removed  through  the  application  of  the  Inverse  of  the  Kalman 
28 

Filter  algorithm.  Program  modifications  are  underway  so  that  the  length 
of  the  polynomial  can  be  automatically  varied  to  achieve  the  desired 
goodness  of  fit  and  error  distribution  based  on  the  statistical  analyses. 
Currently,  the  local  polynomials  are  the  most  efficient  four  coefficient 
polynomial , 

Z • Cq  + CjX  C2y  + CjXy 

the  simplest  non-trlvlal  case.  Note  that  this  order  node  polynomial. 
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tfhen  combined  with  the  weighting  functions  specified  by  WIT,  allows  for 

the  modeling  of  terrain  forms  with  a wavelength  as  low  as  the  node 

polynomial  spacing — approximately  0.2  kilometers.  This  resolution  has 

/■ 

been  found  adequate  In  all  tests  conducted  thusfar,  but  In  any  case  this 

will  not  be  a limitation  after  the  completion  of  the  previously  mentioned 

modification.  When  these  processes  are  complete,  the  node  polynomial's 

coefficients  can  be  written  out  as  the  database,  or  used  for  grid  generation 

or  automated  contouring.  The  first  step  necessary  In  further  use  of  the 

node  polynomial  approximations  Is  their  combination  with  weighting  functions 

as  specified  by  the  WIT  Interpolation  algorithm.  The  WIT  algorithm  defines 

the  final  model  over  each  square  area  defined  by  the  centroids  of  four 

neighboring  node  polynomial  approximations  as  a single,  different  simple 

function.  Currently,  the  simplest  and  most  efficient  polynomial  ratio 

29 

weighting  functions  are  used.  These  weighting  functions  have  a quadratic 
polynomial  numerator  and  denominator,  which,  when  combined  with  the 
currently  used  linear  node  polynomials,  results  in  a final  model  that  is 
a polynomial  ratio  with  a cubic  numerator  and  quadratic  denominator.  The 
globally  smooth  and  continuous  model  defined  by  WIT  results  In  a simple 
low  order  polynomial  ratio  In  each  of  the  square  areas  defined  by  the 
centroids  of  the  node  polynomials.  These  areas  have  been  optimally 
determined  as  0.05  square  kilometers  for  the  current  software  configuration. 
Grid  generation  and  analytic  solution  for  the  contour  lines  are  very 
efficient  because  of  the  simplicity  of  the  global  model. 
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The  general  characteristics  and  strengths  of  this  digital  terrain 


format  can  be  summarized  as  follows: 

1.  automatic  model  accuracy  definition  through  the  use  of  approximation 
theorems,  statistical  analyses,  and  variable  length  polynomials. 

2.  the  modeling  algorithm  is  ameanable  to  implementation  on  a parallel 
processor  because  of  the  identical  Independent  processing  steps  for  the 
local  polynomial  subareas. 

3.  automatic  removal  of  large  data  errors  via  surface  slope  filtering 
and  inverse  Kalman  Filter  algorithms. 

4.  efficient  model  format  requiring  minimal  core  storage  because  of 
the  sequential  and  Independent  determination  of  the  local  polynomials. 

5.  low  storage  overhead  since  the  horizontal  position  of  the  local 
polynomials  are  implied  by  their  location  in  the  matrix. 

6.  a smooth  and  continuous  globally  valid  model  through  the  use  of 

WIT. 

7.  efficient  geographic  access  of  the  dstabase  because  access  to  a 
horizontal  position  is  given  by  the  structure  of  the  matrix. 

8.  efficient  grid  generation  due  to  the  local  low  order  polynomial 
definition  of  the  global  model. 


Development  Plan. 


1 
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Having  identified  the  Polynomial  Matrix  Format  as  the  most  promising 
compression  technique  for  archival  use  In  a production  environment,  the 
folloving  steps  or  stages  of  development  were  identified. 

1.  Determine  the  largest  areas  accurately  modelable  by  the  lowest 
and  highest  order  node  polynomials  considered  for  use. 

2.  Determine  the  minimum  precision  of  digital  representation  needed 
for  the  node  polynomial  coefficients. 

3.  Determine  the  optimal  digital  data  compression  technique  for  the 
node  coefficients. 

The  largest  area  accurately  modelable  under  all  terrain  conditions  with 
the  polynomial  range  chosen  must  be  determined  so  that  the  automatic 
accuracy  algorithm  will  work.  The  area  size  must  be  chosen  so  that  the 
most  efficient  four  coefficient  model  Is  utilized  most  of  the  time,  and 
small  enough  so  that  the  highest  order  polynomial  available  can  accurately 
model  the  roughest  areas  encountered.  Knowing  the  minimum  precision 
necessary  for  the  digital  storage  of  the  node  polynomial  coefficients  will 
permit  the  assignment  of  the  minimum  number  of  bits  per  coefficient. 
Determination  of  the  optimal  digital  data  compression  technique  for  the 
node  coefficients  will  provide  for  the  reasonable  minimum  total  digital 
storage  requirements.  The  reasonableness  of  the  approaches  will  be  judged 
based  upon  such  factors  as  compaction  achieved,  data  regeneration  processing 
time  and  power  requirements,  and  the  effects  of  partial  data  losses. 

The  following  sections  report  the  status  of  our  ongoing  development 
effort. 
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tlmum  Modeling  Area  Determination. 


The  need  for  this  study  Is  based  upon  the  following  empirical 

observation.  The  size  of  the  square  area  optimally  modeled  by  a final, 

weighted  polynomial  when  a four  coefficient  node  polynomial  Is  used 

depends  on  the  terrain  roughness. This  result,  coupled  with  the 

fact  that  this  modeling  algorithm  requires  that  the  size  of  the  square 

area  remain  constant,  dictates  that  node  polynomials  of  varying  complexity 

be  used.  Unfortunately,  research  has  also  shown  that  the  overall 

modeling  efficiency  Is  Inversely  related  to  the  length  of  the  node 

polynomial.  However,  this  penalty  can  be  minimized  if  Dr.  Rauhala's 

30 

Array  Algegra  algorithm  Is  used.  For  this  reason,  the  approach  taken 
here  Is  to  limit  our  analyses  to  consideration  o-f  the  blquadrlc  as  the 
largest  node  polynomial.  Other  factors  contributing  to  this  decision 
are  as  follows.  The  current  optimum  square  ground  area  for  the  final 
polynomials  Is  0.05  square  kilometers  (based  on  four  coefficient  node 
polynomials) , over  a 664  square  kilometer  test  area  (the  area  covered 

c 

by  the  1:50,000  Cache,  OK  sheet)  which  contains  all  representative 
terrain  types.  The  optimum  was  empirically  determined  as  the  largest  square 
areas  which  resulted  In  90  percent  of  the  modeled  points  being  within  - 5 
meters  (one-half  the  contour  Interval)  of  the  stereocompilation  data 
(extreme  errors  excluded) . The  contours  produced  from  this  model  are 
too  generalized  In  the  rough  terrain  areas,  fit  the  noise  In  the  flat 
areas,  and  are  good  in  the  rolling  terrain  areas.  Note  that,  although 
90  percent  of  the  elevation  values  are  within  one-half  contour  Interval, 

90  percent  of  the  contours  are  not!  Further,  three  and  four  coefficient 
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node  polynomials  have  been  investigated  and  there  was  good  correspondence 
between  the  nund>er  of  coefficients  used  In  the  node  polynomials  and  the 
size  of  the  square  area  modeled  to  within  a given  accuracy.  The  best  com- 
paction previously  reported,  30  to  1,  was  based  upon  use  of  the  three 
coefficient  node  polynomials.  Vlhen  the  length  of  the  node  polynomials 
was  Increased  from  three  to  four  the  area  covered  to  within  the  same  accuracy 
Increased  so  that  the  node  polynomials  are  160  times  as  snarce  resulting  In  a 40 
to  1 compression  ration,  (12500-4  coefficient  polynomials  per  2 x 10^  elevations). 

When  these  arguments  are  all  combined  they  result  In  a reasonably  good  | 

expectation  of  modeling  square  areas  as  large  as  .14  square  kilometers 
with  a single  final  polynomial.  This  Is  based  on  the  assumption  that 
four  or  less  coefficient  node  polynomials  could  be  accurately  used  In 
the  flat  areas  and  that  nine  coefficient  node  polynomials  would  accurately 
model  the  roughest  terrain  encountered.  Work  has  been  Initiated  on  Imple- 
mentation of  these  strategies  and  results  of  the  use  of  these  hypotheses 
should  be  available  by  the  end  of  this  fiscal  year. 

I 

I 
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Minimum  Precision  Determination 


At  the  start  of  this  study  the  precision  with  which  the  floating 
point  node  coefficients  must  be  stored  to  maintain  model  accuracy  was 
unknown.  The  only  Information  available  was  that  the  terrain  modeling 
software  had  been  successfully  run  on  large  scale  digital  computers 
with  36  and  60  bit  word  sizes.  In  order  to  Investigate  the  effects  of 
varying  coefficient  precision,  four  measures  of  model  accuracy  were 
Identified;  the  mean  of  the  differences,  the  normalized  root  mean  square 
of  the  differences,  the  standard  deviation  of  the  differences  and  the 
contour  line  locations.  The  mean  of  the  differences  Is  the  mean  value 
of  the  model  minus  grid  elevation  values  and  the  root  mean  square  of 
the  differences  Is  the  square  root  of  the  summed  squared  differences 
divided  by  the  number  of  samples.  The  standard  deviation  has  Its  normal 
meaning. 

The  test  area  chosen  is  a well  surveyed  644‘t-  square  kilometer  region  ' 
covered  by  the  1:50000  Cache,  Oklahoma  map  sheet.  As  mentioned  in  the  last 
section,  previous  empirical  testing  had  shown  that  the  current  modeling 
software  achieved  a modeling  accuracy  of  + 5 meters  for  90X  of  the  stereo 
compilation  data  (over  2 million  points).  The  histogram  of  residuals, 
model  values  minus  measured  values,  for  this  case  la  shown  in  Figure  1 and 
the  model  statistics  are  contained  In  Table  2. 

At  present  the  only  precision  test  completed  Is  a simple  tnincatlon 
Integerlzatlon  of  the  node  polynomial  coefficients.  The  histogram  of 
residuals  which  resulted  are  plotted  In  Figure  2 and  the  model  statistics 
are  contained  In  Table  2.  Figure  3 shows  the  histogram  of  the  differences. 
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I The  results  of  this  first  simple  test  Is  suprlslngly  quite  good.  As 

would  be  expected  from  truncation  Integerlzatlon  the  mean  changed  from 

r 

essential  zero  to  essentially  one  half  meter,  this  shift  should  be  almost 
completely  removed  If  rounding  the  numbers  to  produce  Integers  were 
\ used  Instead  of  trtmcatlon.  The  normalized  root  mean  square  residuals 

differ  by  only  .078  meters  and  the  standard  devaltlons  by  only  .038  meters. 

The  contour  line  plots  were  Indlstlngulshably  different. 

These  are  very  critically  Important  results  In  the  consideration  of 
the  feasibility  and  compaction  of  the  polynomial  matrix  method.  All  com- 

r 

; parlsons  of  compaction  to  date  assumed  that  the  polynomial  coefficients 

would  require  no  more  precision  than  that  currently  used  for  the  elevation 
values,  namely  16  bits.  If  the  node  polynomial  coefficients  had  to  be 
stored  In  a floating  point,  e.g.,  mantissa  and  exponent,  format  then  not 

J 

only  might  this  have  resulted  In  greatly  Inbreased  data  compression  and 

I regeneration  times  but  also  might  have  required  more  than  16  bits  per  ' 

I i 

i 

coefficient.  This  could  have  seriously  adversely  effected  the  compression 
, ratio  for  this  technique. 

The  precision  requirements  were  further  Investigated  by  a statistical 

i 

study  of  the  coefficients  generated  over  the  Cache  test  area.  The  means 
and  standard  deviations  obtained  are  summarized  In  Table  3.  The  terrain 
‘I  over  this  test  area  varied  from  320  to  740  meters.  The  mean  of  the  Co  term 

can  be  Interpreted  as  the  average  terrain  elevation  since  It  Is  the  average 

1 i 

of  the  elevation  predicted  by  the  12,500  node  polynomials  at  their  centroids. 

The  C^  and  C2  term's  mean  can  be  Interpreted  as  the  average  slopes,  and 
I the  Cg  term's  mean  can  be.  Interpreted  as  the  average  cross-slopes. 
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or  interaction  of  the  X and  Y slopes,  e.g.,  this  term  can  be  understood 
as  either  the  change  in  X slope  with  movement  along  the  y dxis  or  vice 
versa.  The  importance  of  these  numbers  is  in  their  magnitude.  The 
magnitude  of  the  term  is  determined  by  the  range  of  elevation  values 
in  the  area  modeled,  and  so  the  use  of  16  bits  for  this  term  is  as 
reasonable  as  the  current  use  of  16  bits  for  the  elevation  values  in  the 
stereocompilation  grids.  The  other  three  coefficients  are  bounded  by 
+2^  ■ ± 128  so  that  if  only  their  Integer  value  is  used  el^t  bits  storage 
for  these  is  more  than  adequate.  This  is  an  important  result  because  our 
compaction  ratio  of  40  to  1 was  based  upon  each  coefficient  being  represented 
with  16  bits.  By  lowering  the  required  length  of  3 of  the  4 coefficients 
the  compaction  ratio  achievable  becomes:  12500  polynomials  with  1-16  bit 
coefficient  and  3-8  bit  coefficients  equals  500,000  bits  divided  into 
2 million  grid  elevation  values  at  16  bits  each  gives  a compaction  ratio 
of  64  to  1. 
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Opt<—i  Digital  D«t«  Co«pr«>alon  Techniques 

Am  stated  before  the  data  compression  philosophy  chosen  in  this  study 
is  that  all  significant  data  smoothing  will  be  accomplished  in  the  mathe- 
matical modeling  step  via  application  of  the  rigorous  least  squares  criteria. 

For  this  reason  only  errorless  digital  data  coiiq>resslon  algorithms  will 
be  considered  for  further  compaction  of  the  Integerized  mode  polynomial 
coefficients. 

The  techniques  identified  for  initial  investigation  is  the  differential 
pulse  code  modulation,  (DPCM),  technique.  This  technique  was  chosen  as 
the  first  to  be  investigated  because  it  is  one  of  the  simplist  both  con- 
cepttially  and  in  terms  of  implementation  difficulty.  The  DPCM  technique 
Involves  storage  of  only  the  difference  from  the  last  parameter  value  en- 
countered. For  example,  the  string  of  nuiid>ers: 

10,  12,  13,  15,  lA,  12,  11,  9 

could  be  stored  using  DPCM  as  the  first  value  plus  the  differences  from 
I the  previous  value,  e.g. 

10,  +2,  +1,  +2,  -1,  -2,  -1,  -2 

The  important  thing  to  note  is  that  the  magnitude  of  numbers  encountered 
In  the  second  case  is  significantly  reduced  and  therefore  wo\ild  require 
fewer  bits  to  represent  them  In  the  digital  computer.  In  order  to  study 
the  applicability  of  the  DPCM  technique  to  the  polynomial  matrix  coefficients 
the  statistics  of  the  absolute  value  differences  between  neighboring  poly- 
nomlala  ware  computed.  The  absolute  value  differences  were  used  because  it 
is  the  Mgnitude  of  the  differences  which  is  of  the  greatest  importance  in 
this  ease.  differences  %fere  computerd  along  rows,  columns  and  diagonals 
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and  the  results  are  in  table  4. 

The  ouch  smaller  range  of  values  for  the  Cq  differences  when  compared 
to  its  original  range  Is  very  significant.  The  largest  number  encountered 
Is  within  "H  ~ 128  so  that  8 bits  could  be  saved  by  using  the  DPCM 
I technique  for  the  Cq  coefficient.  The  range  of  values  of  differences  for 

and  are  essentially  equal  and  all  exceed  the  2^  **  128  range 
that  bound  the  values  of  the  coefficients  themselves.  Therefore,  based 
upon  these  numbers,  use  of  the  DPCM  technique  for  the  C^^,  C2  and  C^ 
coefficients  does  not  appear  to  be  beneficial.  These  numbers,  however. 
Include  some  areas  of  dense,  large  magnitude  errors  which  are  negatively 
affecting  the  range,  so  that  further  analyses  will  be  performed. 

If  the  DPCM  technique  is  used  for  only  the  Cq  coefficient  then  It  can 
be  represented  with  only  eight  bits  Instead  of  16,  this  results  In  a new 
compaction  ratio  of  80  to  1,  since: 

12500  node  polynomials  with  4-8  bit  coefficients  each  equals  400,000 
bits  divided  Into  the  2 x 10^  - 16  bit  elevations  gives  the  80  to  1 ratio. 

Other  techniques  are  being  Investigated  tihlch  take  advantage  of  the 
low  mean  and  standard  deviation  absolute  value  differences.  The  fact  that 
these  are  absolute  value  difference  means  coupled  with  their  relatively 
small  magnitude  Indicates  that  the  nund>ers  are  tightly  grouped  at  the  low 
end  of  the  scale.  Histograms  of  the  differences  were  prepared  which  showed 
that  the  magnitude  of  90Z  of  all  the  differences  were  within  twice  the 

i 

standard  deviations. 
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f Conclusions. 

I 

I The  charscterlstlcs  end  uses  of  automatic  stereocompilatlon  equipment 's 

^ terrain  elevation  data  have  a significant  effect  on  the  suitability  of  the 

various  digital  terrain  formats,  DTF.  The  most  Important  requirements 
of  the  DTF  and  Its  generating  algorithms  are  that  It  a)  must  remove  large 
magnitude  errors,  b)  smooth  the  non-negllgable  frequency  noise  by  making 
optimum  use  of  the  redundancy  In  the  data,  c)  allow  for  efficient  geographical 
access  of  the  data,  d)  permit  efficient  grid  generation  for  use  In  the 
applications  and  c)  result  In  good  compression  with  reasonable  regeneration 
penalties. 

The  various  DTF  approaches  have  been  Identified  and  compared  with  the 
result  that  the  Polynomial  Matrix  method  Is  the  most  promising  at  this  time, 
i A plan  for  development  of  production  algorithm's  and  software  has  been 

formulated  and  work  Initiated.  The  preliminary  results  have  proven  that 

I 

cosqtactlon  ratios  of  80  to  1 can  be  guaranteed  and  much  higher  ratios  can 
I be  expected. 

The  major  tasks  remaining  are  the  evaluation  of  the  applicability  of 
standard  errorless  digital  data  compression  techniques  and  an  analysis  of 
the  tradeoffs  between  compaction  achieved  and  the  costs  associated  with  com-  ! 

I 

^ paction  and  regeneration  for  use. 

i , ! 

I I 
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Table  1.  Digital  Terrain  Format 'a  Characteristics 

(*based  upon  projections  from  insufficient  data 
see  text) 


Normalized 

Standard  Root 

Mean  Deviation  Mean  Square 

(meters)  (meters)  (meters) 

60  BIT  FLOATING  POINT  COEFFICIENTS  0.043  3.809  3.8087 

TRUNCIATED  INTEGER  COEFFICIENTS  -0.558  3.8467  3.887 

DIFFERENCES  0.601  0.0377  0.0783 

Number  of  elevations  * 2 x 10® 


Table  2.  Cache  Test  Area  Model 

Characteristics  Comparison 


COEFFICIENT 


MEAN  STANDARD  DEVIATION  MINIMUM  MAXIMUM* 


6.1 

9.2 

0.0008 

122.2 

^2 

5.5 

8.8 

0.0002 

105.9 

S 

5.1 

7.4 

0.0010 

96.9 

Table  3.  Node  Coefficient  Statistics 
' for  Cache  Test  Area  (based  on 

Model  Equation  - C^+Cj^X+C^T+C^XY) 


NOTE;  Statistics  shown  are  absolute  value  statistics  for  Cj^,  C2t  end  C3 


*This  data  contained  areas  of  dense  large  magnitude  errors  which  were 
not  removed  prior  to  processing. 
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COEFFICIENT 

DIFFERENCE  DIRECTION 

MEAN 

(meters) 

ABSOLUTE  VALUE  STATISTICS 

ST^ARD  DEVIATION  MINIMUM 

(meters)  (meters) 

MAXIMUM* 

(meters) 

(a) 

Row 

5.3 

8.0 

0.000736 

122.3 

(b) 

Column 

4.9 

7.9 

0.000005 

99.5 

(c) 

Diagonal 

6.8 

10.1 

0.001285 

117.5 

W) 

Averages 

5.7 

8.7 

0.001169 

113.1 

C 

1 (a) 

Row 

8.5 

12.4 

0.000241 

164.1 

(b) 

Column 

8.0 

13.4 

0.000283 

192.6 

(c) 

Diagonal 

8.5 

12.2 

0.000271 

146.4 

(d) 

Average 

8.3 

12.7 

0.000265 

167.7 

c. 

^ (a) 

Row 

5.9 

11.5 

0.000065 

156.4 

(b) 

Column 

6.4 

10.9 

0.000487 

124.9 

(c) 

Diagonal 

6.6 

11.0 

0.00071Q 

118.6 

(d) 

Average 

6.3 

11.1 

0.000264 

133.3 

(a) 

Row 

7.2 

11.1 

0.000200 

163.6 

(b) 

Column 

7.2 

*10.7 

0.000262 

151.2 

(c) 

Diagonal 

7.2 

10.4 

0.000276 

128.6 

(d) 

Average 

1 

7.2 

10.7 

0.000246 

147.8 

Table  4.  Cache  Test  Data  Coefficient 
Difference  Statistics 


i i 


*thla  data  contained  areas  of  dense  large  ugnltude  errors. 
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